Effect of Dresselhaus spin-orbit coupling on spin dephasing in asymmetric and 
macroscopically symmetric (llO)-grown quantum wells 
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We develop the microscopic theory of electron spin dephasing in (llO)-grown quantum wells 
where the electron scattering time is comparable to or exceeds the period of spin precession in the 
effective magnetic field caused by spin-orbit coupling. Structures with homogeneous and fluctuating 
Rashba field, which triggers the dephasing of electron spins aligned along the growth direction, are 
analyzed. We show that the Dresselhaus field, which is always present in zinc-blende- type quantum 
wells, suppresses the spin dephasing enabling very long spin lifetime of conduction electrons. The 
dependence of the spin lifetime on the electron mobility is found to be nonmonotonic reaching the 
minimum in structures where the scattering time is comparable to the period of spin precession in 
the effective magnetic field. 

PACS numbers: 72.25.Rb, 75.70.Tj, 78.67.De 



I. INTRODUCTION 



The spin dynamics of free electrons in zinc-blende- type 
quantum wells (QWs) is largely determined by the elec- 
tron subband spin-orbit splitting which can be fruitfully 
treated as a wavcvector-dcpendcnt effective magnetic 
field acting upon the spins of individual electrons [l|-|Jl. 
The field occurs due to the lack of space inversion sym- 
metry in a particular QW and can originate from bulk, 
structure, or interface inversion asymmetry, or their joint 
action. The field direction and magnitude depend on the 
QW crystallographic orientation and design as well as 
external gate voltage |5l-[l2l| , which opens a way towards 
the spin manipulation and spin-based semiconductor de- 
vices. Of special interest are QWs grown on the (110) 
substrate. The effective magnetic field caused by bulk 
inversion asymmetry (Dresselhaus field) in such struc- 
tures points along the growth direction and, therefore, 
does not cause the precession of electron spins oriented 
along the QW normal [l[. This enables the observation 
of very long electron spin lifetimes up to hundreds of 
nanoseconds at low temperatures which are unfeasible 
for two-dimensional structures of other crystallographic 
orientations [l3l-[l7l. 



The electron spin lifetime in (llO)-grown QWs at low 
temperature is likely to be limited by the spin preces- 
sion in the homogeneous Rashba field present in asym- 
metric QWs or unavoidable fluctuating Rashba field that 
emerges due to domain structure formation or inhomoge- 
neous distribution of charged impurities in doped struc- 
tures [l^. The Rashba field lies in the QW plane and 
starts up the spin dephasing. The previous theoret- 
ical studies of electron spin dephasing in (llO)-grown 
QWs with the homogeneous [l^, [20| or fluctuating |2l| - 
[23j Rashba field were limited to the collision-dominated 
regime of spin dynamics, where the scattering time is 
much smaller than the spin precession period in the effec- 
tive magnetic field. Such a consideration is inapplicable 
for QWs with high electron mobility and/or strong spin- 
orbit coupling, which are now technological available. In 



this paper, we develop the microscopic theory of electron 
spin dephasing for arbitrary ratio between the scatter- 
ing time and the spin precession period. We show that 
the interplay of Rashba and Dresselhaus fields leads to 
dynamic coupling of the in-plane and out-of-plane spin 
components. The dependence of the spin lifetime on the 
electron mobility is nonmonotonic and reaches the mini- 
mum in QW structures where the scattering time is com- 
parable to the period of spin precession in the effective 
magnetic field. In high-mobility QWs, the Dresselhaus 
field efficiently suppresses the spin dephasing allowing 
one to reach very long spin lifetimes. 



II. UNIFORM RASHBA FIELD 

First, we consider asymmetric (llO)-grown QWs with 
a spatially homogeneous Rashba field. The linear in the 
wavcvector spin-orbit splitting of the electron subband 
in such QWs (Cg point group) is generally described by 
three linearly independent parameters ai, a2, and /3 |24l | 
(see also Refs. [2^, [2g). The corresponding Hamiltonian 
has the form 



Hso = aiUxky - a2(Jykx + P(Jzkx 



(1) 



where Oj {j 



x,y,z) are the Pauli matrices, kj are the 
components of the electron wavevector. x \\ [110] and 
y II [001] are the in-plane axes, and z \\ [110] is the QW 
normal. The parameters ai and a2 are nonzero in QWs 
with structure inversion asymmetry only, while (3 requires 
bulk inversion asymmetry and describes the Dresselhaus 
splitting in (llO)-oriented QWs. The difference between 
«! and a2 is usually small because it is caused by an 
interference effect between cubic structure of the crys- 
tall lattice and QW structure inversion asymmetry |24| . 
Therefore, we assume in what follows that ai = a2 = a. 
The Hamiltonian ([T]) is similar to the Zeeman term with 
the effective Larmor frequency 



Kp Kp kp 



(2) 





FIG. 1. (a) Dependence of the Larnior frequency flk corre- 
sponding to the effective magnetic field on the wavevector k 
at the Fermi circle, (b) The (xz) plane containing fl^ for 
all the wavevectors k. Precession of electron spins initially 
aligned along the QW normal z with the frequency flk leads 
to the emergence of an average spin component along y. 



where Qji ~ 2akF/h and 51 £> = 2l3kp /h are the frequen- 
eies corresponding to Rashba and maximal Dresselhaus 
splitting at the Fermi level, and kp is the Fermi wavevec- 
tor. Following the experiments [3, [a, [TEI, [131 ^^ consider 
the degenerate electron gas and study the spin dynamics 
at the Fermi level. 

Figure [IJi shows the distribution of the frequency flf, 
on the wavevector at the Fermi level. The distribution is 
strongly anisotropic; the frequency ilk contains both the 
in-plane and out-of-plane components. However, one can 
notice that the frequencies Q,^ corresponding to various 
wavevectors all lie in a certain plane {xz), see Fig[T|D. The 
{xz) plane is obtained from {xz) by the frame rotation 
around the a; axis with the angle 6 = arctan(rifl'/r2i3). 



A. Time dependence 
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FIG. 2. Time dependence of the components of the total 
electron spin S{t) oriented along the QW normal 2: at i = 0. 
Curves are plotted for r2_R = Q.d and different parameters 
Q,Dr characterizing the electron mobility. Dotted curves cor- 
respond to ballistic transport, r = co, and are plotted after 
Eq. (O. Dashed curves present the results of the collision- 
dominated approximation, Q,dt <^ 1, and are calculated after 
Eq. (Uni). 



is obtained by summation Sf,{t) over the wavevectors at 
the Fermi level, which yields 



s,{t) 



E 

j' 






S,.{Q) 



The fact that the plane of the Jlfc vectors does not 
contain the growth direction z nor is perpendicular to z is 
a feature of low-symmetry QWs and leads to the coupling 
of the in-plane and out-of-plane spin components. This 
is illustrated in Fig. [TJd which shows the precession of 
electron spins initially oriented along the QW normal. 
In the ballistic regime, i.e., the absence of scattering, the 
time dependence of electron spin at each point k on the 
Fermi circle is given by 



Sk{t) 



^k ■ s{0) 



^l 



n. 



nk-s{o) ^ ■ 

s(0) ^ Sifc 



^k X s(0) . 

sm{ilkt) , 

\ik 



cos(r2fci) 



(3) 



where s{0) is the initial spin. The time evolution of the 
components of the total electron spin 



S{t)=J2sk{t) 



(4) 



E 

i' 



^l %' 



^ 'k,j ^ ^h,j' 



^l 



cos(r!fcf))S'jv(0), (5) 



where the angle brackets denote averaging over the di- 
rection of k. The first term on the right-hand side of 
Eq. ([5]) presents the averaged projection of individual 
electron spins onto the effective field direction and is 
time- independent. In (IIO) QWs, this contribution is 
present even for the spin polarized along the QW normal, 
which is in contrast to (OOl)-grown structures where the 
effective field lies in QW plane and the time-independent 
term vanishes for S{0) \\ z. The second term on the 
right-hand side of Eq. ([5]) describes spin oscillations. 

Time dependence of the out-of-plane and in-plane com- 
ponents of the total spin oriented along the QW normal 
at i = is shown by dotted curves in Fig.[2|i and[2]3. One 
can see that, in asymmetrical (IlO)-grown QWs, the Sz 
and Sy projections are dynamically coupled: the dephas- 
ing of spins polarized along the QW normal leads to the 
emergence of the in-plane spin component. Both Sz{t) 
and Sy{t) contain oscillations which slowly decay in time. 



Analytical expressions for Sz{t) and Sy{t) at i > 1/^k 
can be derived by the stationary phase method, which 
gives 
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where fly^ = \/^Ji + ^jj ■ The oscillation frequencies 
are determined by the smallest (Hr) and the largest 
(ils) magnitudes of the effective Larmor frequency at 
the Fermi circle. At t — >■ cx), the oscillations vanish and 
the spin polarization reaches a stationary value. 

The scattering of electrons by static defects or phonons 
in real structures leads to the isotropization of spin dis- 
tribution in fc-space which drastically suppresses the spin 
oscillations. In the presence of scattering, the spin dy- 
namics can be described by the kinetic equation 

dsk , _ Sfc - (sfc) 

dt 



SkX^k=g 



(7) 



where g is the spin generation rate at the Fermi level, e.g., 
due to resonant optical pumping with circularly polarized 
light, and t is the scattering time. We assume that g is 
independent of the wavevector direction. The last term 
on the right-hand side of Eq. ([7]) presents a simple form 
of the collision integral. Effects of scattering anisotropy 
on spin dephasing were studied in Ref. [27|. 

Sohd curves in Figs. [5^ and[2b show the time evolution 
of the out-of-plane and in-plane projections of the total 
electron spin S{t) oriented at i = along the QW nor- 
mal. The curves are calculated numerically after Eq. ([7]) 
for different scattering times r. As expected, the decrease 
in the scattering time results in a suppression of spin os- 
cillations. Surprisingly, the effect of electron scattering 
on spin lifetime is more complicated. Comparison of dif- 
ferent curves in Fig. [2^ and [2Jd shows that the shortest 
spin lifetime T '^ r is achieved at a certain mobility cor- 
responding to T '^ 1/^fc. At both higher, r > 1/Ofc, and 
lower, T < l/r2fc, mobility the spin lifetime increases as 
T (XT and T ex 1/{^\t), respectively. 

In the collision-dominated regime, ^rt, Qdt ^ 1, one 
can analytically solve Eq. ([7]) and obtain the master equa- 
tion for the total electron spin [1| 



dt 



^cdS , 



(8) 



where G = X^fcS ^-^^^ ^cd is the D'yakonov-Percl' spin- 
relaxation-rate tensor, (Ted) jj' = t \^k^jj' ~ ^k,j ^ 
For (llO)-grown QWs, the tensor Fed has the form 



■cd 





'^l + ^l 
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i £ p ~r ^ ' r) 


nRfiD 





Oflfii) 


mj. 



(9) 



Solution of Eq. ([8]) with G = and the initial condition 
S'(O) II z has the biexponential form 

c u\ - ^flCxp(-7,j,t) + n% exp(-7gt) 

Sy{t) = ^^^^ ^'"P^-^^'l,; exp(-7.t)] ^^(Q^ (^Q^ 

where jy = (2^^ -|- f2|))r/2 and 75 = nj^T/2 are eigen 
values of the T tensor. The functions (jlOp are plotted in 
Figs. [2ji and [JJd by dashed curves and demonstrate the 
perfect agreement with the results of numerical calcula- 
tion for the collision-dominated regime. 



B. Continuous- w^ave excitation 

Besides experiments with time resolution, spin dephas- 
ing is also widely studied at continuous- wave (cw) pump- 
ing. In the regime linear in the pump intensity, the gen- 
eral relation between the steady-state spin S and the spin 
generation rate G has the form 



S = TG, 



(11) 



where T can be referred to as the spin lifetime tensor. To 
calculate T we solve Eq. (O for the stationary case fol- 
lowing the procedure described in Ref. [2g]. We express 
the spin distribution function Sk via its average value 



Sk 



Cfc + T Ofc X Cfc + -r^ r?fc(Ofc ■ Ck) 
1 + ^1^' 



(12) 



where ^^ = (sk) +gT, sum up Eq. (|12p over the wavevec- 
tor and obtain the linear equation set 



y^ / ^fc 'T Sj]' 



i 'ft; j' i ^k,j' 



1 + ^1^' 



iS,, + G,,T) = G,. (13) 



Solution of the equation set p3)) for the effective Larmor 
frequency ^ yields 



(14) 



v/1 + ri^rVl + f^l^' - 1 



nlji + nlr^) 



ni 














Q, 



ni 



nj,r' 



^ 
^l 



Equation ([T4|) is one of the principal results of the paper. 
It enables one to calculate the steady-state spin polariza- 
tion for arbitrary ratio between the scattering time r and 
the precession periods l/fl/j, I/VLd and shows that the 
Dressclhaus field considerably affects the electron spin 
dephasing in (llO)-grown QWs. 
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FIG. 3. Dependence of the steady-state spin components Sz 
and Sy on the Dresselhaus field for cw spin pumping along the 
QW normal z. Solid curves are plotted for different fiflr and 
normalized by Gz/{Q,\t) which determines Sz at fio = 0. 
Dashed curves correspond to the collision-dominated regime, 
Q.R,Q.D < 1/r. 



Figures [3^ and |3}d show the dependence of the S^ and 
Sy projections of the steady-state spin, caused by exci- 
tation along the QW normal z, on the Dresselhaus field. 
The curves are calculated after Eqs. pT|) and P^ for dif- 
ferent ri/jT. Despite the fact that the Dresselhaus field 
points along the QW normal, it efficiently suppresses the 
spin dephasing increasing Sz as well as leads to the emer- 
gence of Sy. The Sz projection (Fig.|3^) is determined by 
the Tzz component of the spin lifetime tensor. The lat- 
ter equals to l/(J7^r) at Vto = and increases with the 
Dresselhaus field as (Vlr)/^R)\J^ + ^^t"^ at VIjjt ^ 1. 
Such an increase in the spin lifetime is caused by sta- 
bilizing action of the Dresselhaus field on the normal 
spin component. The in-plane spin projection (Fig. Wp) 
is determined by the off-diagonal component Tyz which 
nonmonotonously depends on the Dresselhaus field. 
This component linearly increases with the Dresselhaus 
field as — r2_D(l -I- ri^T^)/(fi|jT) at small VLd^ reaches 
a peak value of -t{1 ^n^T'^fl'^ /[n^.r'^ y/2 + V^] at 
W^ = VLr^{1 -I- ri|,T2)(2 + r2|,r2), and then saturates to 
-(I + n\T'^)/VLB, at VLdt > 1. 

The dependence of the steady-state projections 5*2 
and Sy on the scattering time r is shown Fig. 21 In 
low mobility structures, the spin polarization decreases 
with the increase in r, in accordance with the the- 
ory of D'yakonov-Perel' spin relaxation in the collision- 
dominated regime. In particular, the spin lifetime tensor 
components determining 5*2 and Sy in this regime arc 
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FIG. 4. Dependence of the steady-state spin components Sz 
and Sy on the scattering time r for cw spin pumping along the 
QW normal z. Curves correspond to different ratio Q.b./^d- 



given by Tzz « {n\ + ny)/[nl^{Sl\ + nl/2)T] and Tyz « 
— 0£i/[rifl(ri|j,-|-rj|)/2)r]. At further increase in the scat- 
tering time, the spin polarization reaches a minimum and 
then rises linearly with r, see Fig. [4] Such a behavior 
is caused by the transition from the collision-dominated 
to oscillatory regime of spin dynamics, sec Fig. [2] In 
the oscillatory regime, the Tzz and Tyz components are 
given by Vi]jT /{nRy/fl^~\~Vl^) and -^dt / yJVl^^T^, 
respectively. We note that the nonmonotonic dependence 
Sz{t) is a feature of low-symmetry QWs and absent, 
e.g., in (001) QWs where Sz oc l/r in both collision- 
dominated and oscillatory regimes. 



III. FLUCTUATING RASHBA FIELD 

Now, we analyze the spin dephasing in macroscop- 
ically symmetric QWs where the Rashba field arises 
due to QW domain structure or inhomogeneous distri- 
bution of charged impurities in the doping layers [l8|, 
[29l . [SOf and fluctuates in the QW plane around its 
zero average value. In such systems, the Larmor fre- 
quency Jlfc corresponding to the total effective field 
is given by the sum of ^D,k = ^d(0, 0, fc^/fc^?) and 
position-dependent Jl_R,fc(p) ~ VLji{p){ky/kp, ~kx/kp,0) 
with p being the in-plane coordinate. The fiuctuating 
Rashba field is usually weak, therefore, we assume that 

\nRip)\<^\nDli/T. 

During the spin dephasing, electrons travel in the QW 
plane and their spins are affected by the Rashba field 
of different strength and sign. The characteristic time 
of the Rashba field change for an individual electron Tc 



(correlation time) is determined by the domain size and 
electron mobility. Depending on the ratio between r^, and 
other relevant times, different scenarios of spin dephas- 
ing are realized. Below, we consider them for the initial 
condition of spin polarization along the QW normal. 

(i) Large domains, Tc » 1/{Vi\t) > \/{VI\,t), t. Spin 
dephasing occurs in all domains independently. To cal- 
culate the time evolution of electron spin in each domain, 
we solve kinetic Eq. (O with 9 = and the initial con- 
dition Sfc(O) II z. At J7fl;/57£) ^ 1, the z axis is nearly 
parallel to the z axis, see Fig. [Ijb). As the result, the in- 
plane component of the total electron spin in the domain 
is small. The out-of-plane component exhibits exponen- 
tial decay 



^.(t)=5z(0)exp(-7jt) 
with the dephasing time 



Iz 



(15) 



(16) 



by 



The net electron spin averaged over all domains is given 



Sz{t)^ / s,{t)w{nR)dni 



(17) 



where 'w{Q.fi) is the distribution function of the Rashba 
field in the QW plane. Different dephasing rates 75 in 
domains lead to a non-exponential decay of the net elec- 
tron spin. In particular, at large times {t 3> I/t? but 
still t ^ Tc), the time dependence of the net spin has the 
form 



SS)=w{0) 



7:{l + ^l + nlT^) 



Tt 



1/2 



(18) 



The presence of the Dresselhaus spin-orbit coupling leads 
to a slowdown of spin dephasing in all domains and, 
hence, to an increase in the net spin polarization. 

(ii) l/(ri|.T) > Tc > l/(fi|,T), T. The process of spin 
dephasing occurs in different domains at different rates 
7j, however, electrons frequently move from one domain 
to another. Therefore, the time decay of the net electron 
spin is monoexponential 



^,(i)=5,(0)cxp(-7?t) (19) 

and described by the domain-averaged rate 



Iz 



Qlip)T 



i + ^i + nir^ 



(20) 



(iii) SmaU domains, 1 / (nj^r) , 1 / {QJ^t) '> Tc :> t. In 
this regime of spin dephasing, electrons move from one 
domain to another more frequently than the 75 and 7^ 
rates. The in-plane component of the net electron spin 



does not emerge due to the frequent change in the Rashba 
field sign. The out-of-plane component exponentially de- 
cays at the rate 



r. 



nl{p)T 



(21) 



which is independent of the Dresselhaus field strength. 
The dephasing rate (PT|) is twice as high as the rate ([20)1 
even at ^dt < 1- 

(iv) l/(i7|.r) ^ T ^ Tc- The correlation time of the 
Rashba field is shorter than the scattering time, and spin 
dephasing occurs at ballistic electron transport. Such 
regime is realized in high-mobility QW structures where 
the fiuctuating Rashba field stems from the electric field 
of charged impurities in remote doping layers (iSl |23|. 
The spin dynamics of an electron with a certain wavevec- 
tor k is described by 



— — + Sfc X [ilD,k 

at 



^rMPo +vt)] = 0. 



(22) 



where po in the initial electron coordinate, v = vpik/k), 
and vf is the Fermi velocity. Equation ([^^ can be solved 
by changing the spin coordinate frame to the one rotating 
along the z axis with the angular frequency fln.k- The 
procedure is equivalent to the substitution s^ = Rk Sfc 
and ^R^k ~ Rk ^B..k with Rk being the rotation matrix. 



Rk — 



cos[il]:)t{kx / kp)] — sm[ilDt{kx / kp)] 

sm[ilDt{kx/kF)] cos[ilDt{kx/kF)] 

1 



This yields dsk/dt + Sk x QR^k{po + vt) ~ 0. Taking into 
account that the Rashba field lies in the QW plane and 
CIr{p)tc ^ 1 we obtain the dephasing rate of the z spin 
component for electrons with the wavevector k 

00 

r,4fc) = / nii^kipo + vt) ■ fiR,k{po) dt (23) 



^r{P{) + vt) Qr{po) coslflntikx/kp)] dt 



where the overline denotes averaging over the initial elec- 
tron coordinate. 

The spatial correlation function of the Rashba field 
caused by Coulomb impurities randomly distributed in 
the (5-doping layers positioned at the distance Zd below 
and above the QW has the form (see Ref. [18|) 



nR{p)nR{p') = 



»fl(p)(2^rf) 



[(p - p')2 + (2zd)2]3/2 

Combining Eqs. (P^ and (P^ one obtains 



(24) 



T,,{k) = niip) nDTl{kx/kF)Ki[nDTc{kx/kp)] , (25) 
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FIG. 5. Dephasing rate of the out-of-plane spin component 
in QWs witli tlie fluctuating Rashba field in various regimes. 



where Tr is the correlation time of the Rashba field 



/ n,R{po + vt) ^r{po) dt 





Kipo) 



2zd 



(26) 



and Ki is the modified Bessel function of t he secon d kind. 
At floTc ^ 1, the rate rzz{k) is given by 17^(p)tc. 

The dephasing rate of the net electron spin in the 
structure can be obtained by averaging T^zik) over the 
wavevector direction, which yields 



(r,4fe)) = nl{p)T,^MOKiiO - hi^KoiO] , (27) 

where /^ is the modified Bessel function of the first kin d 
and £, = Q.£,Tc/2. The dephasing rate is given by ri|j.(p)Tc 
(see Ref. [l3|) and il|,(p)/r2D in the cases of moderate, 
^dTc ^ 1, and strong, VIdTc ^ 1, Dresselhaus field, 
respectively. 

The theoretical results on spin dephasing rate in quan- 
tum wells with the fluctuating Rashba field are summa- 



rized in Fig. [5] The strong Dresselhaus field increases the 
electron spin lifetime in all regimes of spin dymanics. 



IV. SUMMARY 

We have developed the microscopic theory of electron 
spin dephasing for asymmetric and macroscopically sym- 
metric (llO)-grown quantum wells where the Rashba field 
is homogeneous or fluctuates in the quantum well plane, 
respectively. Both low-mobility and high-mobility struc- 
tures are analyzed. It is shown that the regime and 
rate of spin dephasing are determined by the strength 
of Rashba and Dresselahaus fields, the latter being al- 
ways present in zinc-blende-type structures, as well as 
the electron mobility. The Dresselhaus field points along 
the growth direction in (110) quantum wells and increases 
the lifetime of the out-of-plane spin component. 

In asymmetric structures, the interplay of the homoge- 
neous Rashba and Dresselhaus fields leads to a coupling 
of the in-plane and out-of-plane spin components. The 
spin lifetime of the out-of-plane component nonmonoton- 
ically depends on the electron mobility. It linearly de- 
creases with the scattering time r in low-mobility struc- 
tures, reaches the minimum at r comparable to the pe- 
riod of spin precession in the effective magnetic field, and 
then linearly increases with r in high-mobility structures. 

In quantum wells with the fiuctuating Rashba field, the 
scenario of spin dephasing is determined by the scale of 
the field fiuctuations compared to other relevant scales. 
We show that, for all regimes of spin dynamics, the strong 
Dresselhaus field slows down the dephasing rate of the 
out-of-plane spin component. 
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